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Abstract 

The initial value problem is considered in the present paper for bipolar quan- 
tum hydrodynamic model for semiconductors (QHD) in M'^. We prove that 
the unique strong solution exists globally in time and tends to the asymptot- 
ical state with an algebraic rate as t ^ +00. And, we show that the global 
solution of linearized bipolar QHD system decays in time at an algebraic de- 
cay rate from both above and below. This means in general, we can not get 
exponential time-decay rate for bipolar QHD system, which is different from 
the case of unipolar QHD model (where global solutions tend to the equihb- 
rium state at an exponential time-decay rate) and is mainly caused by the 
nonlinear coupling and cancelation between two carriers. Moreover, it is also 
shown that the nonlinear dispersion does not affect the long time asymptotic 
behavior, which by product gives rise to the algebraic time-decay rate of the 
solution of the bipolar hydrodynamical model in the semiclassical limit. 

Key words: Quantum hydrodynamics; Algebraic decay rate. 



1 Introduction 

The quantum hydrodynamic(QHD) model for semiconductors is derived and studied re- 
cently in the modelings and simulations of semiconductor devices, where the effects of 
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quantum mechanics arises. The basic observation concerning the quantum hydro- 
dynamics is that the energy density consists of one additional new quantum 
correction term of the order 0{s) introduced first by Wigner | |29j in 1932, and 
that the stress tensor contains also an additional quantum correction part |[2J, (3j 
related to the quantum Bohm potential (or internal self-potential) [i4j 

with observable p > the density, m the mass, and e the Planck constant. The 

quantum potential Q is responsible for producing the quantum behavior. Such possible 
relation was also implied in the original idea initialized by Madelung [25] in 1927 to 
derive quantum fluid-type equations in terms of Madelung's transformation applied to 
wave function of Schrodinger equation of pure state. Recently, the moment method is 
employed to derive quantum hydrodynamic equations for semiconductor device at nano- 
size based on the Wigner-Boltzmann (or quantum Liouville) equation, see in [26] for 
details. For derivation about quantum hydrodynamical equations and related quantum 
models, one can refer to [6l [71 [15] and the reference therein. 

In this paper, we consider the Cauchy problem of the bipolar quantum hydrody- 
namic(QHD) model for semiconductors in x [0, +oo) which reads 

dtPi + V ■ (piUi) = 0, (1.2) 

^tip^Ui) + V ■ ip^u, ® Ui) + VPiipi) = q^p^E + ^p,V(— ^) - (1.3) 

2 Jpi Ti 



AV ■ E = Pa- pb- C{x), V X E = 0, E{x) ^ 0, l^l +00, (1.4) 
with the initial conditions 

{pi, Ui){x, 0) = (pi„, UiJ(x), (1.5) 

where the index i = a,b and qa = 1, qb = —1- The variables pa > 0, pb > and Ua,Ub 
and E are the particle densities, velocities and electric field, respectively. We can define 
the usual momentum Jq, Jb as Ja = PaUa, Jb = PbUb- Pa{-) and Pb{-) are the pressure- 
density functions. The parameters e > 0, Ta = u = t > 0, and A > are the scaled 
Planck constant, momentum relaxation time, and Debye length respectively. C = C{x) is 
the doping profile function. When it holds (pt, PbUb) = (0, 0) formally, the above model 
reduces the unipolar quantum hydrodynamical model. 

Recently, many mathematical efforts are made on the study of the QHD model for 
semiconductors on both the steady state solutions and the evolutional (time-dependent) 
solutions. The investigation on unipolar QHD model are well-understood up to now. 
The steady state solutions of unipolar QHD model are studied in [5l [ini [IE], [30] in one- 
dimensional or multi-dimensional bounded domain for different boundary conditions, and 
the steady state solution of the unipolar viscous quantum hydrodynamical system is 
investigated in |j8j. For the one-dimensional time- dependent case, the short time existence 
of solutions of unipolar model [TT] and the global existence theory with the exponential 
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stability of stationary state in whole space [TU [TTl [12] are established. For the multi- 
dimensional case, the local existence of solutions is obtained for irrotational fluid [19], 
and the local and global existence theory and exponential stability of equilibrium state 
analysis are also investigated for irrotational fluid on spatial periodic domain [22]. The 
corresponding existence theory for time-dependent solution for general rotational fluid 
is usually difficult and is obtained very recently in [13], where the exponential decay to 
the stationary state obtained therein is made. Moreover, the asymptotical small scaling 
analysis including the relaxation time limit, small Debye length limit and the semi classical 
limit for the global solutions are studied in [TSl EI] respectively. 

However, the results for bipolar QHD model are quite fewer compared with those 
obtained for unipolar QHD model. So far, only the steady state solutions are studied 
partially in [201 |28] for bounded and unbounded domain, and the semiclassical limit 
and relaxation limit of the global-in-time solutions are investigated in [31], where the 
global existence of time- dependent solution is also proven, but without the deriving the 
large time behavior. The main difficulty in dealing with the bipolar QHD model is the 
coupling and interaction between the two carriers, which may cause some cancelation, 
and it is not clear that the equilibrium state to the bipolar QHD is still exponential stable 
or not for small perturbation. 

In this paper, we study the time-decay rate of global solutions to the Cauchy problem 
for the bipolar QHD ([L2!)-([L5D in We shall show that the solution to the IVP for 
bipolar QHD tends to the equilibrium state at an algebraic decay rate. This property is 
different from the unipolar QHD model and is caused by the interaction and nonlinear 
coupling of the two carriers which make the convergence of solution to the equilibrium 
state slower. 

We have the following main result. 

Theorem 1.1 Assume C{x) = c* with c* a positive constant, and p* > 0, pi > are 
constants satisfying pi - pi - c* = 0. Assume Pa, Pb e and P'SpD-, PUpD > O.Let 
the initial data satisfy {p^^ — p*,Uig) G H^iM.^) x H^iM.^), i = a,b, with Aq := \\{pio — 
Pi 5 ""io) ll-H''^(R^)xW5(]j3). Then, there exists Ai > such that if Aq < Ai, the unique 
solution {pi,Ui,E) of the IVP <\1.2^ - (\1.5f\ with pi > exists globally in time 
and satisfies for i = a,b that 

{p^-p*) e C\0,T-H'^-^\R^)), Ui e C''{0,T-n'^-^\R^)), 

E eC\0,T;n^'^''(R^)), (1.6) 

for k = 0,1,2. 

Moreover, the solution {pi,Ui,E) tends to the equilibrium state (p*,0,0) at an 
algebraic time- decay rate 

(1 + t)'\\D'{p, - p*)f + (1 + tnsD^p, - p*)f < cAo, < ^ < 5, (1.7) 
(1 + t)''\\D''{ui, Ji)f + (1 + t)''\\D''Ef + (1 + t)^\\D^Ef < cAo, 1 < A; < 5, (1.8) 
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where the coefficient c> is independent of e, and 7i^(M'^) denotes the space 
that {/ G L^{R^),Df e H''~^{R^)}, k>l. f denotes the k-times spatial deriva- 
tive of f. 

Remark 1.2 By ( |i. 7[ )- pTgj) and Nirenberg's inequality for three dimensional case 

11 i 1 

||«||l-(R3) < c\\D\\\l,^^,^\\u\\le(^^,^ < C||D\|||2(K3)||^«|I12(K3) (1-9) 

we can get the optimal L°° time- decay rate of the solution 

\\{p,-p*,u„E)\\L^^u^)<c{l + t)--^. (1.10) 

This time- decay rate is the same order as the heat equation in three dimension. In fact, 
when taking relaxation limit for bipolar QHD, we can get the bipolar quantum Drift- 
Diffusion (QDD) equation below. For this bipolar QDD model, we can show 
that the global solution of initial value problem tends to the equilibrium state 
with the same rate as heat equation ^24j. 

dtPi + V ■ [qiPiE - VPiipi) + ^p.V(^^)] = 0, (1.11) 

^ \ Pi 

A^V ■ E = Pa- pb- C{x), V X E = 0, E{x) 0, |x| ^ +oo. (1.12) 

Unlike the unipolar quantum hydro dynamical model [HI [22l [121 US], in Theo- 
rem [HT] we can not get the exponential convergence to the asymptotical equilibrium state 
for bipolar quantum model for the whole space case due to the coupling and cancelation 
interaction between two carriers. In fact, by the original equations fll.2p - fll.4l) . we can get 
the linearized system around the equilibrium state for the variables 

{Wa,Ja,Wh,Jb,E) = {pa - p*^, PaUa, Pb - pi, PbUb, E) 

that 

' Wat + V-Ja = 

Jat + P'aipD^Wa - ^VAl^, + J„ _ p*^ = Q 

< Wbt + y-Jb = Q (1.13) 

Jbt + P'bipD^Wb - f VAWfc + J, + plE = Q 
^ V ■ E = Wa-Wb V X E = 0, E ^0 as l^l ^ oo 

with initial data given by 

{Wa, Ja, Wb, Jb)ix, 0) = {WaO, JaO, WbO, Jbo){x) (1.14) 

where we have let r = 1, A = 1 for simplicity. From the Poisson equation 01.131) ^; for the 
electric potential E we can represent E by 

E = VA-\Wa-Wb). (1.15) 
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Assume that the initial data (11.141) satisfies 

Jao, Jbo G H'{R'), Wao, Wko G H%R^) n L\R^), (1.16) 

so that the initial electric field Eq obtained from Poisson equation (ll.lSp at initial time 
has the regularity 

Eo = VA-\Wao - Who) G H'{M.'). (1.17) 



Remark 1.3 The norm \\D''E\\ of E with integer k > can be obtained by 
Lemma \2.1\ through the Poisson equation and the norm \\E\\i2 is from the 
Riesz's potential theory in that \\E\\ < c\\(Wa — Wb)\\ e with a positive con- 
stant c. 

For simplicity, we just consider the IVP (I1.13I) - (I1.14I) for following case 

^ = 1, P: = 2, = 1, c* = 1, P:(2) = Piil) = 1, (1.18) 

since the method used in section 4 to prove theorem 11.41 about the time-decay rate of 
solutions to IVP f ll.13p - fll.14p can be applied to general case instead of fll.l8p . 

We have the algebraic time-decay rate of global solution to IVP problem (11.130 - 01.140 
for the case fll.l8p below. 

Theorem 1.4 Suppose that i\1.16^ - i\l.lS\) hold. Assume that the Fourier transformation 
{Wao, Wbo) of initial density satisfy for some constants mo > 0, r > that 

inf J(#ao + 2iy,o)(OI >^o, (1.19) 

and the initial perturbation of momentum satisfies 

V-(J,o + 2Jw) = 0. (1.20) 
Then, the unique global solution to { \1.13\i -( 1.14'^ exists and satisfies 



Wa,Wb e Ci[0,+oo),H%W')), Ja,Jb e Ci[0,+oo),H\R')), 
E eCi[0,+oo),H\W')), 

and 

ci(l + t)-t-t < \\{d',Wa,d^,Wbm\\L2^u:^) < C2(l + t)'^ < A; < 6, (1.21) 
ci(l + t)-M < \\{d',Ja,d^,Jbm\\LHm < C2(l + t)"l < A; < 5 (1.22) 

fori = a,b. The positive constants Ci,C2 depend on thq, ||f/o||//6x//5, and \\{Wao,Wbo)\\L^ ■ 
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Remark 1.5 The theorem \1.4\ shows that for above linearized bipolar QHD, the density 
and momentum have only algebraic time- decay rate from both above and below. This fact 
means that in general one can only expect an algebraic time-decay rate for the original IVP 
problem for nonlinear bipolar QHD 1.4), since the nonlinear bipolar QHD system 

can be viewed as a small perturbation of the corresponding linearized system. 

As one can see that all the estimates (11. 71) - (11. 8 1) and (ll.lOp hold uniformly with respect 
to the Planck constant e, thus we can apply the theorem established in [31] to pass into 
the semiclassical limit e — 0+ in ( 11.2p -( fL5|) . and obtain the algebraic time decay rate of 
the following limiting solution (which is the solution of the limiting equation- the classical 
bipolar hydrodynamical model) as £ — 0+ below 

dM+W ■{p,u,) = Q, (1.23) 

Oil- 

dMu,) + V ■ (p,u, ® Ui) + VPiipi) = QiPiE - (1.24) 

A^V ■ E = Pa- pb- C{x), V X E = 0, E{x) 0, \x\ +oo. (1.25) 

We have the following result about the decay rate of the corresponding solution of 
bipolar HD model as an application of Theorem 11.11 in the process of semiclassical limit. 

Theorem 1.6 Under the assumptions of Theorem there exists {pi,Ui,E), i = a,b, 
such that as e 0+, the solution {p^,u^,E'^) of IVP tends to {pi,Ui,E) 

strongly 

Pl^p, m c{Q,T-clnHf-% u\^u, m c{Q,T-cl nn1;,% 

E^' ^E m C{Q,T-Ctnnl:), sG(0,|). 

And where {pi,Ui,E) is the solution of the bipolar HD model U.2!^) - D.25\) with initial 
data ^1.5\) . Moreover, it also holds 

II (p, - pIu„ E) |Uoo(k3) < c(l + t)-t (1.26) 

as t ^ +00. 



The rest part of the paper is arranged as follows. After some preliminary given in 
section 2, we shall prove Theorem 1.1 and Theorem II. 61 in the section 3, and we will prove 
Theorem 1.3 in Section 4. 



2 Some preliminary 

Notations 

C and c always denote the generic positive constants. L^(R^) is the space of square 
integral functions on M? with the norm || ■ || or || • ||l2(]r3). H^{E?) with integer k > 1 
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denotes the usual Sobolev space of function / satisfying dlf G L^(M'^)(0 < i < k) with 
norm 

0<|a|<fc 

here and after a G N'^,D°' = d^^d^^d^^^ for |a| = si + S2 + S3, Especially || ■ ||o = || • ||. 
Let i3 be a Banach space, ^^([0, t]; B) denotes the space of i3- valued k-times continuously 

differentiable functions on [0,t]. We can extend the above norm to the vector- valued 

3 

function u = {ui,U2, ,U3) with = ^ and 

r=l 



D'ur= [ {Y,Y.^D-Urf)dx, 



k 

and \\u\\k = ||M||Hfc(R3) = \\D'u\\, \\f\\L^([Q,T];B) = sup \\f{t)\\i3. We also use the space 

i=0 0<t<T 

H*^(M3) = {/ g L%R'^),Df e H^-\M.^)},k > 1. Sometimes we use || (.,.,...) ||//fc(K3) or 
||(., ., ...)||fe to denote the norm of the space H''{R^) x H''{R^) x ■ • ■ x H''{R^) and the 
n''(R^) as well. 

Lemma 2.1 Let f G H^(R^), s > |. There is a unique solution of the divergence equation 
V ■ u = f, VxM = 0, u{x) — > 0, |x| — ^ +00. 

satisfying 

\\u\\l(^(M.3) < C||/||2,2(K3), ||D-U||^S(]R3) < C\\f\\Hs(R3y 

We will also use the Moser type calculus lemmas. 
Lemma 2.2 Let f,ge H'{R^)f]L°^{R^), then it holds 

< C\\g\\L^ ■ + CII/IUo. . \\D^g\\ 

WD'^ifg) - fD^gW < C\\g\\L^ ■ + ■ {{D^'^^-'gW 

for a G A^^, 1 < \a\ < s, s > is an integer. 

Lemma 2.3 Let f G H''^{R^) with s > be an integer and function F{p) smooth enough 
and F(0) = then F{f){x) G H'{R^) and 



3 The proof of Theorem 11.11 and Theorem 11.6 



Note that the local and global existence of the solution in Theorem 1.1 can be referred 
to [31] , we only focus on the convergence rate of the solution to the corresponding steady 
state. 
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3.1 The reformulation of original problem 



Our idea is to obtain the uniform estimates of the local solution, and we need to refor- 
mulate the original problem into a convenient form. Take A = 1, r = 1 and use {.)t to 
denote dt{.) for convenience. First, by equations fll.2p -( |T73f) we can get the equations for 
i^i = \fPi = ^) as in [31] 



0, 



(3.1) 



with the initial value 



ipit{x,0) := ^pii{x) 



By equation (11.31) with the fact [ui ■ V)ui = |V(|ttip) — Mj x (V x Ui), taking curl 
of the two sides of the equation (II. 3p we get for 0j = V x as 

(Pu + 0i + (Ui ■ V)(P^ + <PiV ■ Ui - {(f)i ■ V)Ui - u,{V ■ 0) = 0. (3.2) 

ipi — a/p* then the 



Here we have V ■ = 0. Introducing new variables Wi 
system for {wa, Wb, (l)a, 4>b, E) is 

Watt + Wat + + -^{Wa + VP^) V • E - P'^{pl)AWa = fal, 

wm + wu + — ^ - -{wb + ^)V • E - Pl{pl)Awb = hu 

(pat + (pa = fa2, 
4>bt + 4>b = fb2, 

V-E = wl-wl + 2y^^Wa-2y^,Wb, V X E = 0, 
with the initial conditions given by 

Wi{x, 0) := Wio{x) = ipio - a/p*, 0j(a;, 0) := (pia{x) = V x Ui^ix), 
Wit{x, 0) := Wii{x) = (-Mi„ ■ Vwj,, - ^{wig + ^/p*)V ■ WiJ 



(3.3) 

(3.4) 

(3.5) 
(3.6) 
(3.7) 



(3.8) 
(3.9) 



and where 

fil ■= fiiix,t) 



-wt. 



Wi + Jp* 



qNw.E + {P[{{w,, + ^,f) - P:{p*))Aw, 



*^2^ 



\Vw^ 

Wi + J7i 
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/i2 := fi2ix,t) =((0, ■ V)Mi - (Ui ■ V)0, - ■ tli), (3.11) 

for i = a,b. We will also use the relation between V ■ Ui and Vwi, Wu from (1.2) 

2wit + 2ui ■ Vwi + {wi + ^/p*)V ■Ui = 0. (3.12) 

3.2 The a-priori estimates 

Assume that the classical solutions Wi,Ui, E satisfy a-priorily 

6t = max {E(l + + E(l + t)'\\D%f + ^(1 + t)'+'\\D'wu\\' 

k=0 k=l k=0 

3 

|2 



+ 5^(1 + ty+'WD'uuf + (1 + tfWD'wuf 

k=l 

5 2 

+ 5^(1 + t)'=||D'=£;f + 5^(1 + t)3+'=||DWf } < 1. (3.13) 

k=l k=0 

It follows for the sufficiently small 5t the positivity of density -^j {i = a, b) as 

^ 3 ^ 

By Nirenberg's inequality for three-dimensional case from fl3.13p . we have 

3 2 
^(1 + tf^^WD'^w^W^ + 5^(1 + tf^^'WD'^w^tW^ + (1 + t)i«;.«||ioo < cbT (3.14) 

fc=0 fc=0 
3 1 3 

5^(1 + t)^'+iD'=^i,iii.. + 5^(1 + t)'^+3||D'=^i,,iii^ + 5^(1 + t)'=+iD'^E|ii^ < c^T. 

fc=0 fc=0 A;=0 

(3.15) 

With the help of the a-priori assumptions (13.131) we establish the following a-priori 
estimates 

Lemma 3.1 For the short time solution {wi, Ui, E) it holds for t G [0, T] that 
5 3 



^(1 + t)''\\D''w,f + (1 + t)5||£DVf + 5^(1 + tf+^WD'wuf 

k=0 k=0 

2 

+ (1 + tnD'wuf + Y^il + tf+'WD'wutf < cAo, (3.16) 



A:=0 
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J2{^+t)''\\D'u.f + J2i^+ty+'\\D'uur<cAo, (3.17) 

k=l k=l 

J](l + t)''\\D'Ef + r ^(1 + s)'-'\\D''Efds < cAo, (3.18) 

k=l k=l 

f {^(1 + sf~'\\D'w,f + ^(1 + s)''+'\\D''wur}ds < cAo, (3.19) 

k=l k=0 

f\j2il + s^D^f + ^(1 + s)'+'\\D''uur}ds < cAo, 

"^0 k=l k=l 



(3.20) 



provided 6t is small enough, where the Aq is defined in Theorem 1.1. 



Proof: Step 1 (the basic estimates). Multiplying equation (13.31) by {wa + '2.Wat), and fl3.4l) 
by (iWfe + '^Wbt), integrating by parts the resulted equations over M^, omitting without 
confusion, summing the resulted two equalities and noticing the fact from Poisson equation 
( K7\i that 

j {(^iwa + v^)V ■ E){wa + 2wat) " {^{wb + v^)V ■ E){wb + 2wbt)}dx 
[ I V ■ E\'^dx +-I I V ■ E\^dx -- I V{wl - wl) ■ Edx 



1 d 
lit 



we can get 



d 
di 



+ ^l{\Aw,\'' +\Awb\'') + ■ E\'}dx 



+ \\V-E\'}dx 



= \j iyu - ^l) ■ Edx 

+ j{fal{x,t){Wa + 2Wat) + fbl{x,t){Wa + 2Wbt)}dx. (3.21) 

By assumptions (13.13^ . using Sobolev imbedding theorem and Holder's 
inequality, Young's inequality and integration by parts, we can estimate the 
right-hand side terms of ( 13.211) as follows 

J WiVwi ■ Edx <l|u'^l|L3||Vu;,|U2||^||i6 
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<c(||«;,|U2 + \\Vw,\\L2)\\Vwi\\L2 ■ WEWls 

<c6Ti\\Vw^f + \\V-Ef), (3.22) 



and 



d f 

-'jtj + " p/(ft-)]iv«).rrfi + *ii(v«„,«,.,)ir , (3.23) 

/,„ . V.,,(2u,,)rfx = - / V . u^^dx < *ll«...f . (3.24) 

j UiV{ui-Vwi) ■2witdx < -— j {ui-Vwifdx + c5T\\{Vwi,Wit)\\'^, (3.25) 

where we have used the fact H-Diijp < c(|| V ■ -Ujp + ||V x UiW^), and ||Vwjp, HifitH^ 
to estimate V ■ Ui through equation fl3.12p . The other terms in the right-hand side of 
f l3.2ip can also be estimated easily by integration by parts, Holder's inequality, Young's 
inequality and the Lemma 2.2 and Lemma 2.3, together with fl3.22l) - fl3.25p we can have 
from (13^11) that 

d f uP' vlP' 

- j {Wl + WaWat + Y + Wl + W^Wu + y + P'M)\VWa\'' + Pfe' (p^ I V^^fe | ' 

+ / iUt + <t) + P'aipDY^^a? + n'(P:)|V^.P + ^(|A^.,p + \Aw,\') + i|V ■ E\^]dx 

< c5T\\{VWa,VWb,Wat,Wbt,V ■ E,0a,0b)f . (3.26) 

Taking inner product between (13.50 and 20^, and between (13. 6p and 20;,, integrating 
over R'^, we obtain 

^ y (I0J2 + |0^|2)^^ + 2 j (|0,|2 + l^feH^x = j {fa2 ■ 20a + fm " 20J(ix. (3.27) 
A simple analysis to the right-hand side of ( 13.270 together with ( 13.130 gives 
^ j{\(t)a? + \M'')dx + 2 y"(|0ar + |06r)^ix<c5T||(0a,0fe,V«;a,Vti;6,u;at,w;6t)||'. (3.28) 
Integrating of the summation of (13.261) and (13.281) over [0, t] and using 

Pl{p*) > 0, ^{x' + y')<x' + xy + ^< 2ix' + y'), 
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we obtain 

II (Wa, Wb) II? + II (eD'wa, sD^Wb) f + II {Wat, Wbt, 0a, <Pb, DE) f 

+ / {\\{Vwa,Vwb,eD'^Wa,eD'^Wb,WauWbtAaAb)\? + \\DEf]ds 
Jo 

< cAo- (3.29) 

Making summation between the integral / { (I3.3P x 2(1 + t)wat + (13. 4p x 2(1 + 
t)wbt}dx and (13.271) x (1 + t), we can have after a comphcated but straightforward 
computation that 

+ t)l {wl + wl + P^ipDlVw^l' + HipDlVwbl' + j{\Awa\' + \Awb\') 

+ + ^bf) - n^PlW-^h? + {Ua ■ Waf + {Ub ■ Wbf}dx} 

+ 2{l + t)\\{Wat,WbtAaAbW 
<{l + t) II {WaU Wbt, (l)a, <Pb) ir + c5t\\ (VWa, VWb, sD^Wa, sD^Wb, 0a, 06, DE) f (3.30) 

where we have used the a-priori time-decay rate assumptions ( 13.131) . Holder's 
inequality, Young's inequality to estimate the right-hand side terms as follows 

+ / {l+t){{fal{x,t){Wa + 2Wat) + fbl{x,t){Wb + 2Wbt)}dx 



+ j {l+t){fa2-2(j)a + fb2-2(j)b}dx 

<{l + t) II {Wat, Wbt, 0a, 06) f + cStW {VWa, VWb, sD^Wa, sD^Wb, 0a, 06, DE) f. 

The integrating of (I3.30p over [0,t] together with the help of (I3.29P gives rise to 



{1 + t)\\{VWa,VWb,eD Wa,eD Wb,Wat,Wbt,4>a,(pb,DE)\\' 
+ / {1 + s){\\{Wat,Wbt,4>a,(pb)\\'^ds 







< cAq. (3.31) 

The combination of (13.291) and (13.311) shows the basic estimates in Lemma 3.1 as 

\\wif + {l + t)\\Dwif + (1 + t)\\DEf + {l+t)\\Duif < cAo, (3.32) 

/ \\{Vwi,DE)fds+ [ {l + s){\\Duif + \\wuf)ds<cAo, (3.33) 
Jo Jo 
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Step 2 (the higher order estimates). Next, we will do the higher order estimates. To this 
end, set Wi := V^Wi, (pi := -D"0i, E := D^Eli = a,b. 1 < \a\ < 4). Differentiating 
equations fl3.3p - fl3.7p with respect to x, we get the equations for Wi,(f)i,E that 



2 

wut + + ^A^^y, - Pl{p*)/\w, + |(«;, + v/p:)V ■ E 



=D"/a(x,t) - D^qiw, + Vp*)V ■ i?) + |(«;. + v/p*)V ■ E, (3.34) 

0i, + 0i = DVi2, (3.35) 
V ■ ^ = D'^iwl - wl + 2^«;, - 2^,w,), (3.36) 

= l,gb = -1. 

Similarly to deriving the previous basic estimates, combining the following 
integrals together 

/ / J2^(^2Mi=ay<{l + s)\D''wa + 2D''wat) + (^3Mi=b^{^ + s)\D''wb+2D"wbt)}dxd^^ 

Jo J 

J J J2{(BB:=a ■ 2(1 + S)'D>„ + (BBi=b ■ 2(1 + s)'D>,}dxrfs 



and 



^ J {mMi=a X 2(1 + sp+'D^Wat + mM^=b X 2(1 + s) l^l+^D^w;, Jdxds 

^ y {mB^=a ■ 2(1 + S)I"I+1D"0, + ^m^=b " 2(1 + Sp+' D"<P,}dxds 

for |a| = A; with k = 1,2,3,4 respectively, we can get after a straightforward 
computation that 

(1 + t)"^' II {D'^^'Wa, D'^+'W,, sD'^+^Wa, sD'^+^W,, D^W^u D^Wu, i^Va, ^ Vb, ^'+'^) 

+ Al + s)'^ II {D^^'wa. D'+'w,, sD'+'w,, sD'+'w,, D'+'E) \\'ds 
Jo 

+ [ {l + s)'^'\\{D'wat,D''wM,D''(Pa,D'(f>,)fds<cAo. (3.37) 
Jo 

By (I3.37P we can get part of decay rates in Lemma 3.1 that 

{l + t)''\\D''w^\\^ + {l + tf\\eD^w4^ <cAo, 0<k<5, (3.38) 

{l + tf\\D''uif + {l + tf\\D''Ef<cAo, l<k<5, (3.39) 

+ sf-^UD'wu D'E)f + (1 + s)'\\D%f}ds < cAo, 1 < A; < 5, (3.40) 

Jo 
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and 

{l+ty+''\\D''wuf+[{l + sy+''\\D''wufds<cAo, 1 < k < 4. (3.41) 

Jo 

The higher order estimate (1 + t)^||D^i?|p < cAq can be obtained by Poisson 
equation ( 13.71) and the Lemma 2.1. 

To complete the proof we still need to do the decay rate of {wa, Wb,Ua, u^) about higher 
order derivatives on time t. Set Wi = D'^wu, (pi = D'^cpu, E = V^Et (0 < |a| < 2), then 
we get the equations for Wi.cpi.E 

2 

Witt + Wit + ^A^w), + |(w, + ^)V • E - P'i{p*)/\Wi 

= D"{fa{x,t))t - D^i^iw, + Vp:)V ■ E)t + ^{w, + ■ E, (3.42) 

0,, + 0, = D"(/,2(x,t))i, (3.43) 
V ■ ^ = D^iwl - wl + 2^^wa - 2^,w,)t, (3.44) 

with i = a,b, Qa = 1, = -1. 

Based on the results derived in fl3.38l) - fl3.41l) we can get from fl3.42l) - fl3.44l) the more 
faster time-decay rate for Wi,(j)i, E as before. Summing the integrals 

nt ^2+|a| 

Yl {(M)«=a(l + sYiD'^Wat + 2D''watt) + (^Mi=b{^ + s)\D''w,t + 2D''w,tt)}dxds 

Jo J 

/ / {^M^=a ■ 2(1 + s)'D^(t>at + ^M^=h " 2(1 + s)' ct>,t}dxds 

Jo J 

and 

y* j {^m^=a X 2(1 + Sp^'D'^Watt + ^m^=b X 2(1 + Sp+^D'^Wttt 

+ <^Mi=a ■ 2(1 + sp+^D'^^at + <^Mi=b ■ 2(1 + sp+''D''(f)bt}dxds 

for a with |a| = 0,1,2 respectively which together with the help of the results fl3.38p - 
fl3.4ip gives us finally 



{l + tp^D'wu\\'+ / {l + sy+'\\D'wurds<cAo, 

Jo 

(1 + tp-'\\D'wut\\' + [\l + sp'WD'wutfds < cAo, 

Jo 

(1 + tp^D%tf + [\l + sp'\\D%tfds < cAo, 

Jo 

{l + tp^D'^Etf + [ {l + sy+''\\D''Etfds<cAo, 
Jo 



0<k<3, 


(3.45) 


0<k<2, 


(3.46) 


0<k<2, 


(3.47) 


l<k<3. 


(3.48) 
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Note that \\Dut\\^ < c{\\V ■Utf) + \\V xutf with the help of I^M)-i^M) and the relation 
of V ■ Mj and Vwi, Wu through the equation (I3.12p . we have 

{l + ty+''\\D''uuf+[{l + sy+''\\D''uitfds<cAo, l<k<3. (3.49) 

Jo 

Then, Lemma 3.1 follows from (13:4511 - (IMjl and (ICTjl - dCTll and (l332l) - ( 13331) . □ 



3.3 The proof of main results 
Proof of Theorem 11.11 and Theorem II. 61 - 

From Lemma 3.1, we know that the sufficiently small Aq makes us be able to extend the 
solution to the global one by continuity argument and the estimates (I3.16p - (l3.20p hold 
for any t > especially that 

(l + t)'=||D'=Wif + (l + t)^||5DVf <cAo, 0<A;<5, (3.50) 

{l + t)^-'^\\D^Witf + il + tf\\D^wuf<cAo, 0<k<3, (3.51) 

{l + t)''\\D''uif + {l+t)''\\D''E\\^ <cAo, l<k<5, (3.52) 

(1 + tf+''\\D''uitf < cAo, 1 < A; < 3. (3.53) 

The coefRcient c is independent of the Planck constant e and time t. As 
Pi = (wi + y/ffY we can get the conclusion of the Theorem 11.11 that 

(1 + tf\\D\p, - p*)ir + (1 + t?\\eD%p, - p*)f < cAo, < A; < 5, (3.54) 

{l + tf\\D''ui\\^ + {l + tf\\D''E\\^<cAo, l<k<5. (3.55) 

Thus, the proof of Theorem ll.ll is completed. From (13. 54^ — ( 13.551) . using Niren- 
berg's inequality we have 



||(a-P*,m„^)|Uoo(k3) <c(l + t)-4. (3.56) 

Let us turn to the proof of the Theorem 11.61 Since all above a-priori estimates es- 
tablished for the solutions given in Theorem 11.11 hold uniformly with respect to Planck 
constant £. Denote the solution by {pl,ul,E^) and it follows that(see[31J) there is a 
solution denoted by (pa, Ua, Pb,Ub, E) such that 

pl ^ A in C(0, T; Cl n Hl-J)- u\ in C(0, T; Cl H Hf"/); 

E-^E in C{0,T- Ct n Hf~/), s e (0, i), 

for any T > 0,i = a,b. One can easily verify that (p^, Ua, Pb,Ub, E) is the global-in time 
solution of the bipolar hydrodynamic model (ll.23p - (ll.25p . What's more, we have the 
estimate by (13.560 that 

-p*,tii,^)|Uoo(K3) <c(l + t)-i (3.57) 
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4 Algebraic decay rate for linearized system 

In this section, we will prove Theorem ll.4[ Namely, we shall show that for linearized 
bipolar QHD system, the density and momentum converge to its asymptotical state at 
an algebraic decay rate from both above and below. This implies that in general we can 
only get an algebraic time-decay rate for bipolar QHD. This is caused by the interactions 
between two carriers. Since the nonlinear bipolar QHD system is a small perturbation 
of the corresponding linearized system, one can only expect the similar results for the 
original problem. 

By fll.lSp . the equation (11.131) for U = (Wa, Ja, Wb, Jb) can be rewritten as 



W„t + V ■ 4 = 

Jat + - VAWa + Ja - 2V A'^Wa - Wb) = 

Wbt + V-Jb = 

Jbt + VWb - VAWb + Jb + VA-\Wa -Wb) = 
with initial data given by 

U{x,0) = Uo{x) =: iWao,Jao,Wbo,Jbo)ix). 
Let us write the solution of the linear problem fl4.ip - fl4.2l) formally 

U = e^*[/n 



(4.1) 



(4.2) 



(4.3) 



where U will be the inverse of its Fourier transformation U = {Wa, Ja,Wb, Jb) whose 
equation can be derived by taking Fourier transform with respect to x on (14.11) as 



where the Matrix 



A 



Ut = AU 

U{^,0) = {WaO,JaO,WbO,Jbo) 

( -i^* \ 

-i^* 
V \id^ -1^62 -h j 



(4.4) 



with 



6l = l + |eP + -q^, C?l = -q7' ^2 = l + |el'+^, ^2 = ^, /3 = diag(l,l,l) 



and the notation i is the imaginary unit. Here J a = {Ja'\ Ja'\ Ja'''), Jb = {Jl^' , J'b', J^ )■ 
We solve the O.D.Es (14. 4p straightforward by linear O.D.Es theory and get its solution 
denoted by 

(j = e^'Uo (4.5) 



j{l) j{2) ^(3)^ 
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where U = {Wa, Ja, Wb, Jb) with 

Wa{^, t) ^^WaoiFi + 2F2 + er + e+ + 2(62 + €+)] 
+ ^WbolFi -F2 + e^ + et- (e^ + e+)] 

- ^ilo ■ 0(^1 + 2F2) - '^{Jw ■ 0(2i^i - 2F2), (4.6) 

WhiC, t) =lWbom + F2 + 2(er + e+) + {e^ + e+)] 
+ lWao[Fi - F2 + er + e+ - (63 + e+)] 

- ^(J,o • 0(2i^i + ^2) - ^(Jao • 0(^1 - ^2), (4.7) 



and for k— 1, 2, 3 



j{k) . 3 

ifne,t) =^(ier-eDe-*- J^(E^'^S)^-* 



i=i 



(^ ■ 4o)[2F2 + Fi - 2(62" + e+) - (e+ + ef)] 



6 

+ J^i^ ■ Jbo)[F2 - Fi + + ei 



+ 



m 



- 2{Wao - Wbo)^F2 - ^WaoU^ + \C\''){2F2 + F,) 
-JWboUl + \^f){F^-F2), (4.8) 



—^{i ■ Jm)[F2 + 2Fi - (62- + 62+) - 2(e+ + e^)] 



where 
with 



m 

+ ^(^ • ■/ao)[i^2 - Fi + er + e+ - 62- - e+] 

+ (W^ao - Wb,)'^F2 - \Wu,U^ + |er)(i^2 + 2Fi) 

-^^aOefc(l + len(i^l-i^2) (4.9) 

e- = e-t(i-^^), e+ = e-i(i+^^), 62" = e't^i-^^), 62+ = e-i(^+^^) (4.10) 

h = ^i-Ami+m, /2 = vi-4(3 + ie|2(i+ie|2)), (4.11) 
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and 

= , F, = . (4.12) 

Note here that we have = V A-^{Wao-Wbo) G L'^{M.^) which imphes (Wao-Wbo)^ e 

L^(M^) in fl4.8p - fl4.9p . This means the existence of the inverse transformation of If and 
thus the global solvability of U for flilD-flO). 



Proof of Theorem\1.4\ We first focus on the estimates of the lower bound in fll.2ip - 
fll.22p . The idea is to analyze the Fourier transformation of U due to the Plancherel 
theorem. In view of fl4.10p - fl4.12p we should give some properties of the terms contained 
in Wa,Ja,Wt, Jb given by dUD-gH). 

We have the following estimates 



+ 1 + \et\ + 


|i^2| + 1^1^21 <ce-=*, 


for 




(4.13) 




ieiVil<cK^*, 


for 




(4.14) 






for 




(4.15) 



where and below c > is a generic positive constant. The estimates ( I4.13h is 
gained by a direct computation. The estimates ( 14. 141) . (14. 151) can be obtained 
as follows. 

It holds for > ^ that 

i-4i^r(i+ier)<o, /i = vi - 4KP(i + m = wmw +m-i = w 

and 



By 



i^ii = i J 1 = 12'''^ — iin — ^' '''''^ l(OI<l, 

we know 

\Fi\ < |e-3*. 

As for leiVil, it holds for ^ < < ^ that 

ierii^il<4e-^*- 

When ^ < leP, we can directly compute 
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By the fact that 1 - a/1 - 4s(l + s) < 2(^2 + l)s for < s < we can obtain 

( KWf easily since = e-t^- Vi-4l€P(i+l5P)) > g-^I^P*. 

With the help of fl4.13l) - fl4.15l) we can turn to calculate the time-decay rates of density 
and momentum Wa, Wb, Ja, Jb, and we take Wa, Ja for simplicity. Set 



Wa = Ti + Ri, Ji^^ = rf') + A: = 1, 2, 3 



with 



Ti = liWao + 2Wto){F, + er) (4.16) 
Ri = Wa-Ti (the rest terms) (4.17) 
= -^iWao + 2^0)6(1 + |er)i^i (4.18) 

^(fc) _ j{k) _ J.{k) ^^^^ ^^^^ ^^^^g^ ^^_^g^ 

By fl4l3l) -f l4:T5l ) we know 

^ Jr3 Jr3 

|2jf „/i I j.\„~ct 



2y|5i2<v^ 6 

>]: I \\{Wao + 2mo)e~A^di-c{l + t)e 

2 J|f|2< v^-^ 



> / ce-^'^l^l - c(l + Oe-'^* 

'|C|2<min{^|=^,r2} 



2 

3 



> c(l + t)-5 -c(l + t)e-'* (4.20) 

where we have used the assumption in Theorem 11.41 that \{Wao + 2Wfco)| > ?t^o > in 
S(0, r) and the fact Fi > for < We also used (l^rW'ao, l^rW^w, l^l'^ao, |^|'^) e 

L2(/?3) for the integers < n < 6, < / < 5 and ^i?M to get I^J^c?^ < ce~^* with 
the help of (14.131) . The above c > denotes the generic positive constant depending on 
the norm of initial data and mo and not necessarily be the same. 

The combination of Plancherel theorem and inequality fl4.20p implies for t ^ 1 that 

\\Wa{.M = Wa{.M>Cl{l+tr^ (4.21) 

with ci some positive number. Similarly, with the help of (I4.13l) - fl4.15l) . we have 
||iW(.,t)f >^ / \iCkT,\^dC- [ HkRil'd^ 

[ \^{Wao + 2Wto){F, + e^)\^d^-c{l + t)e-''' 

2 ./|^|2< b 



2 
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> c [ |efene-2^l«l'*|de-c(l + t)e-^* 

i|g|2<mm{^|^,r2} 

> c(l + t)-5 -c(l + t)e-"*. (4.22) 
It follows from (141221) that 

\\d,,Wa{.,t)f = \\i^kWa{.,t)f > Ci(l (4.23) 

for t 3> 1. Repeating the similar procedure as above, we can estimate the higher order 
term ||i'"'i^"^^2^i^3^M/a|p (|a| < 6), which together with the Plancherel theorem leads to 
the algebraic time-decay rate for Wa from below 

\\diWai.,t)\\L2^RS)>cil + t)-'2-l, 0</<6. (4.24) 

Again, we can repeat the similar argument as above to establish the corresponding al- 
gebraic time-decay rate for J^. In fact, by f l4.13p -f HrT5l) we have after a direct computation 
that 

\\Ji'\:t)r>l I iTf r^e- / \Rf\'di 

>\ I |Tf rrfe-c(l+t)e-* 

>c[ l^kFil^d^- €{1+1)6-"', fc = l,2,3. (4.25) 

i|^|2<mm{^^,r2} 

Note that < Ji < 1 for < I^P < |(^), we have 

li^il = \- 1 = ^|e-5(i-^i)(l - e-"^)\ > ce-5(i-^i) (4.26) 

h \h\ 

for t > 1 and \(\^ < |(^). By KWf and the fact e-K^-^i) > ce'^l^l'* for \^\^ < |(^), 
we finally obtain from (14.261) that 

>ce-'=l«^ for |e|2<l(^). (4.27) 

Set rf = min{r2, |(^)} and let t > 1. By Km . g^ZD, we have 

\\j!,'\.,t)r>c [ iefc|V2^i«i'*rfe-c(i + t)e-^* 

>c{l + t)~l-c{l + t)e-"\ k = 1,2,3. (4.28) 

This gives rise to the time-decay rate of Ja = {J!i'\ J'P) for t ^ 1 that 

||j.(.,t)|| = ||j.(.,t)||>ci(l + t)-t. (4.29) 
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The higher order estimates of Ja can be estabhshed in the similar argument as obtain- 
ing (14.221) for Wa and finally we can have for t ^ 1 that 

\\DlJai.,t)\\ > c(l + / = 1,2,3,4,5. (4.30) 

The above estimates are valid for Wt, Jh^\k = 1, 2, 3) due to the symmetry between Wa 
and Wb, Ja and J^. Thus the proof of the lower bound estimates in Theorem ll.4l is finished. 

Note that the time-decay rate from above of solutions in fll.2ip and (11.221) can be 
obtained in the same framework of Fourier transformation to establish the lower bound 
of decay rate. Also, it can be obtained by energy methods used in Section 3, we omit the 
details. □ 
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